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^vq I PACS 47 . 55 . D- - Drops and bubbles 

PACS 68 . 08 . Be - Wetting 

f II ' Abstract. - Liquid droplets impacting a superhydrophobic surface decorated with micro-scale 

O . posts often bounce off the surface. However, by decreasing the impact velocity droplets may land 

^ ■ on the surface in a fakir state, and by increasing it posts may impale droplets that are then 

stuck on the surface. We use a two-phase lattice-Boltzmann model to simulate droplet impact on 

e"""^ , superhydrophobic surfaces, and show that it may result in a fakir state also for reasonable high 

' impact velocities. This happens more easily if the surface is made more hydrophobic or the post 

height is increased, thereby making the impaled state energetically less favourable. 
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^ \ Introduction. — Recently, there has been rapidly growing interest in designing ar- 

tificial surfaces with extreme hydrophobic properties. These efforts have been inspired in 

y^ I particular by biological superhydrophobic surfaces such as plant leaves or insect wings [1,2]. 

■<!:;;j- ' These surfaces have a hydrophobic coating and they typically have micron-scale roughness 

which has been found to further enhance the hydrophobicity of the surface. Similar designs 

f~^ ' have been utilized in order to manufacture artificial, strongly water-repellent surfaces (see, 

00 ; e.g., refs. [3,4]). 

The fact that roughness has a major effect on the hydrophobicity of a surface has long 
been recognized. As a droplet is deposited on a rough hydrophobic surface, it can be found 
in two different states. In the first one, known as the impaled or Wenzel state [5], the 
droplet follows the surface topography and no air is trapped beneath it. The second state 
C^ ' is known as the fakir or Cassie-Baxter state [6] . In this state the droplet sits on top of the 

roughnesses, and some air remains trapped in the hollows and grooves under the droplet. 
For the fakir droplets, a high contact angle is found. Also the contact angle of a droplet 
in the Wenzel state is higher than that observed for one on a smooth surface, although not 
as high as in the fakir state [7]. Usually, a high contact angle is related to a low contact 
angle hysteresis, which is the most important property of a superhydrophobic surface, even 
though this relation is not generally true [3,8,9]. Therefore, in order to obtain strongly 
water-repellent behaviour, surface designs leading to fakir droplets are preferred. 

The fakir state is often metastable and therefore an impalement transition to the Wen- 
zel state may occur as a result of some external disturbance [10-12]. As the fakir state 
is typically the preferred one, it would be important to understand the robustness of this 
state, and the mechanisms leading to collapse to the Wenzel state, which may dramatically 
change the hydrophobicity. This issue has recently attracted a lot of interest, and both ex- 
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perimcntal [11, 13, 14] and numerical [14-16] techniques have been utilized to elucidate the 
transition between different superhydrophobic states. These studies have been conducted on 
microstructured surfaces decorated with regular patterns of posts. The impalement transi- 
tion can be triggered by several mechanisms, including, among others, external pressure [10], 
defects on the surface [11], or hydrophilic contaminants deposited on the surface from the 
air or from the droplet itself [12]. 

One possibility to study impalement transitions is to consider droplets impacting mi- 
crostructured hydrophobic surfaces [17, 18]. In this case the dynamic 'pressure' due to the 
impact velocity of the droplet is the mechanism triggering the impalement transition. Re- 
cently, Bartolo and co-workers [17] reported of three distinct regimes, in which the qualitative 
behaviour of an impacting droplet is different depending on the impact velocity. First, if 
the impact velocity is small enough, droplet lands on the posts and after some oscillations 
stays on the surface as a fakir droplet. Second, for intermediate impact velocities, droplet 
bounces off the surface. Third, for high impact velocities, sticky droplets are observed. In 
this regime, the posts on the surface impale the droplet, and liquid penetrates the volume 
between the posts. Thus, in the last case, the droplet is found in the Wenzel state. 

In this work, we study droplets impacting on a hydrophobic surface patterned by posts 
with square cross section. We are especially interested in the circumstances under which a 
droplet is impaled by posts, and when a non-bouncing droplet is observed. Our approach is 
numerical and utilizes the latticc-Boltzmann method to simulate the behaviour of impacting 
droplets. In particular, we find the same three regimes as reported by Bartolo et al. in 
ref. [17]. In addition, under certain circumstances, we find a fourth possible regime for 
impact velocities larger than those for bouncing droplets but smaller than the ones for 
sticky droplets. In this state, a non-bouncing droplet that ends up in the fakir state is 
observed. 

Lattice-Boltzmann method and the Shan-Chen model. The present simula- 
tions were done using the latticc-Boltzmann (LB) method. As this method is well established 
and throughly covered in a number of review articles and books (see, e.g., refs. [19-23]), we 
only briefly discuss here, for the sake of completeness, the basics of the method. 

In the LB method, the fluid is described by an ensemble of particles moving along 
links between lattice nodes of a regular lattice. Time and velocities are discretized such 
that during one time step particles can move only to neighbouring lattice nodes. The LB 
fluid properties are determined by single-particle distribution functions fi (r, t) which can be 
interpreted as the probabilities to flnd a particle at lattice node r at time t moving with a 
discrete velocity c^. These distribution functions evolve according to the LB equation, 

/,(r + c„f + 1) - fdv,,t) ^ -- {Mr,t) - /r(r,i)) . (1) 

T 

Here, the right-hand side of the equation describes the collisions among the particles as a 
relaxation process towards a local equilibrium which is a low-Mach-number expansion of the 
Maxwell-Boltzmann distribution. We use the relaxation-time approximation with a single 
characteristic time scale r. 

The macroscopic quantities are obtained as velocity moments of the distribution func- 
tions. In particular, the mass and momentum densities are given by 

P = E/^ (2) 

i 

pu = ^c,f, (3) 

i 

respectively. 

In order to model a two-phase fluid, we use the multiphase model developed by Shan 
and Chen [24] . In this model a mean-fleld interparticle interaction is added to the standard 
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LB equation. This force can be expressed in the form 






F = Gb^(r)}_^t,^(r + c,)c„ (4) 

i 

where ip ~ 1 — exp (— p/po) is an effective mass (po is a reference density) and Gb is a 
parameter that adjusts the strength of the interaction. The values of weight factors ti depend 
on the magnitude of the corresponding discrete velocity c^. Here we use the standard D3Q19 
model (three-dimensional lattice with 19 discrete velocities), and the weights for this model 
can be found in ref. [25]. This additional interaction enables simulation of liquid- vapour 
systems with surface tension. The wetting behaviour of the liquid at solid walls is modelled 
with a similar type of force added between the solid and the fluid. Here we set a density 
value Pw to the solid lattice nodes, and the contact angle can be modified through this 
density. 

Results and discussion. — We use the following geometrical setup in the simulations. 
The size of the system is 150 x 150 x 170 lattice nodes, where the last dimension is related to 
the vertical direction. To the bottom of the system we add an array of posts with a square 
cross section of size 3x3. The distance between neighbouring posts is 10 lattice spacings. 
The length of the posts is between 10 and 25 lattice spacings. 

In the middle of the simulation domain a droplet with a diameter of 105 lattice spacings is 
initialized. In the beginning of a simulation, we give the droplet an initial (impact) velocity 
Uo- We choose such values for the simulation parameters that the density ratio of the liquid 
and vapour phases is approximately 30. Contact angle is slightly varied in the simulations, 
and the intrinsic contact angle (i.e., the contact angle of a droplet on a perfectly smooth 
surface) is between 106° and 120°. Notice that all lengths above were given in dimensionless 
lattice units and the results below are given in terms of appropriate dimensionless numbers. 

First we simulated such a case in which the intrinsic contact angle was 113°, and the 
post height was 20 lattice spacings. We varied the initial velocity and measured the contact 
time, i.e., the time that the droplet was in contact with the surface before bouncing off. 
In the non-bouncing cases the contact time obviously approaches infinity. Experimental 
studies [17,26] and scaling arguments [26] suggest that contact time does not depend on the 
impact velocity if this velocity is high enough. However, in their experiments Richard et 
al. found that, when the Weber number (the ratio of fluid inertia to surface tension) was 
We — pL/qR/j <C 1, contact time increased with decreasing impact velocity [26]. With the 
parameters used in the present simulations. We varied between 0.1 and 7. Therefore, the 
contact times shown in fig. 1 behave as expected. For high impact velocities the contact 
time stays essentially constant, but for somewhat lower velocities, when We is decreased 
below unity, it increases for decreasing impact velocity. 

Superficially the contact times shown in fig. 1 lead us to conclude that there are three 
different regimes in the qualitative behaviour of droplet impact. First, when the impact 
velocity is low, the contact time is infinite. Second, there is a regime of an almost constant 
contact time. Third, another regime with infinite contact time is observed at high impact 
velocities. These regimes would correspond to the non-bouncing, bouncing, and sticky 
droplets, respectively, as in ref. [17]. However, a closer inspection of the last regime reveals 
two different types of qualitative behaviour. For very high impact velocities we find the 
sticky behaviour where the posts impale the droplet. But for lower impact velocities when 
the contact time still diverges, the droplet is eventually found in the fakir state. In fig. 2 
we show a series of snapshots from the impact process in this case. It is evident that posts 
impale the droplet, but as the droplet bounces, the interpost volume is drained during the 
process. However, the droplet does not bounce off but remains on the surface and finally 
ends up in the fakir state. We thus find it in a similar state as in the non-bouncing case 
although the two processes are different. We distinguish between these two states by calling 
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Fig. 1: The inverse of dimensionless contact time as a function We ' — t/o\/pi?/7 (impact veloc- 
ity). Time is made dimensionless by scaling it by to — pR^/fi where /i is the dynamic viscosity. 
The intrinsic contact angle is 6 = 113°. 



them the first and second non-bouncing state corresponding to an impact velocity that is 
lower and respectively higher than that of the bouncing droplet. 

Next we considered how the hydrophobicity of the surface affects the occurrence of the 
second non-bouncing state. To this end, we simulated droplet impact on surfaces with 
varying hydrophobicity such that the intrinsic contact angle varied between 106° and 120°. 
The post heigth was fixed to 15 lattice spacings and also the impact velocity was varied so 
as to find the velocity intervals where different impact regimes, and especially the second 
non-bouncing state, are found. The results are shown in fig. 3. We observe that, as the 
hydrophobicity grows (i.e., the contact angle increases), the second non-bouncing state is 
found at higher impact velocities, but also the velocity interval of this state becomes wider. 
Also notice that, for the lowest values of contact angle used in the simulations, we do not 
observe the second non-bouncing state at all regardless of the impact velocity. It is thus 
evident that the second non-bouncing state is found easier if the surface is made more 
hydrophobic. 

In a similar fashion, we studied the effect of post height. The contact angle was kept 
contant, 9 = 112°, and the simulations were performed on surfaces with post heights between 
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Fig. 2; Cross sections of a 3D droplet-impact simulation in the second non-bouncing regime. (1) 
The initial droplet. (2) Droplet impacts the surface and starts to fill the volume between the posts. 
(3) Droplet tries to bounce off but remains stuck at the centre of the contact area. (4) Maximum 
height for the centre of mass during the attempted bouncing. (5) Centre-of-mass motion is again 
downwards but, at the same time, the volume between posts begins to be drained. (6) Droplet in 
its final state. 
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Fig. 3: Impact regimes for varying hydrophobicity (intrinsic contact angle) and impact velocity 
(yWe — Uo \/pR/'y) ■ The regimes are first non-bouncing, bouncing, second non-bouncing (white), 
and sticky regime, from left to right. The simulation results are depicted with stars, filled circles, 
open circles, and crosses, respectively. The post height in all simulations was 15 lattice spacings. 



10 and 25 lattice spacings. The results of these simulations are shown in fig. 4. The 
qualitative observations from these simulations are similar to those reported above. As the 
height of the posts is increased, a wider interval of impact velocities leads to the second 
non-bouncing state. Also, a higher impact velocity is necessary for this state as the post 
height increases, and when the posts are short enough, the second non-bouncing state is 
absent. 

The results described above are related to the metastability of fakir droplets. As already 
discussed in the Introduction, the fakir state is typically metastable, but transition to the 
energetically more favourable impaled state requires external triggering in order to overcome 
the energy barrier separating these two states. Our investigation shows what happens to 
an impacting droplet when the impaled state is made less and less favourable. This is 
exactly what happens when the hydrophobicity of the surface or the length of the posts are 
increased. 

The intrinsic contact angles leading to the second non-bouncing state appears to be 
quite large. However, such values are still experimentally achievable as the largest contact 
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Fig. 4: Impact regimes for varying post height and impact velocity {\/We — UoyJ pR/^). The post 
height is made dimensionless by scaling it by droplet radius R. The intrinsic contact angle in all 
simulations was 112°. The colours and symbols are as in fig. 3. 



angle observed on a smooth surface is presumably about 120° as reported by Nishino and 
co-workers [27]. These authors used trifluoromethyl carbon self-assembled on a surface. 
Another possibility to achieve high contact angles could be to utilize roughness at different 
length scales such that the patterned surface itself would be made rough with some smaller- 
scale roughness. This type of surface design has been proposed by Patankar [28]. 



Summary. — In conclusion, we simulated droplets impacting hydrophobic surfaces 
patterned with regular arrays of posts. Our results were in agreement with those of Bartolo 
et al. [17], i.e., the behaviour of an impacting droplet depends on the impact velocity such 
that a non-bouncing, bouncing or sticky droplet results from the impact as the impact 
velocity increases. We found in addition that, when the impaled state is made energetically 
less favourable, another non-bouncing state is observed between the bouncing and sticky 
state. This state might be useful for gaining better understanding of the energy-barrier 
related mechanisms in droplet impalement and thus for the development of more robust 
superhydrophobic surfaces. 
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